Cavity QED of the graphene cyclotron transition 
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We investigate theoretically the cavity quantum electrodynamics of the cyclotron transition for 
Dirac fermions in graphene. We show that the ultrastrong coupling regime characterized by a 
vacuum Rabi frequency comparable or even larger than the transition frequency can be obtained for 
high enough filling factors of the graphene Landau levels. Important qualitative differences occur 
with respect to the corresponding physics of massive electrons in a semiconductor quantum well. In 
particular, an instability for the ground state analogous to the one occuring in the Dicke-model is 
predicted for increasing value of the electron density. 
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Cavity Quantum Electrodynamics (QED) in the ul- 
trastrong coupling regime is a fascinating topic that 
is attracting considerable interest in condensed matter 
physics, particularly in semiconductor microcavities[ll-@| 
and superconducting circuit QED systems The ul- 

trastrong coupling regime is achieved when the vacuum 
Rabi frequency (quantifying the interaction between one 
cavity photon and one elementary electronic excitation) 
becomes comparable or even larger than the correspond- 
ing electronic transition frequency. In such regime, it is 
possible to manipulate the quantum ground state of the 
cavity system, to modify the decoherence properties of 
the systemjiol and to enhance interesting non-adiabatic 
cavity QED effects 0, [ll|. 

Recently, it was theoretically predicted that by cou- 
phng a cavity photon mode to the cyclotron transition (of 
frequency wo) of a two-dimensional electron gas (2DEG) 
in a semiconductor, it is possible to have a vacuum Rabi 
frequency such that ilo/'^o ^ ^afs v uq^ where afs 
is the fine structure constant, is the filling factor of the 
Landau levels and nqw is the (effective) number of quan- 
tum wells [Til]. Such a predicted scaling has been quanti- 
tatively demonstrated by recent impressive spectroscopy 
experimental results in the THz domain[l3] .The strik- 
ing consequence of this physical behavior is that for high 
filling factors it is possible to have flo/ujQ ^1. In this 
ultrastrong coupling limit for the 2DEG, however, no in- 
stability is expected for the ground state, as a result of 
the role played by the so-called diamagnetic A'^-term[ll| 
(A is the electromagnetic vector potential operator). In 
the case of massive electrons in a semiconductor, the ef- 
fective mass approximation is known to work very well 
in a broad range of conditions and one can generally con- 
sider just the conduction band to describe the physics of 
the two-dimensional electron gas, hence the underlying 
crystal structure of the semiconductor host turns out to 
be unimportant for many physical effects. In the case of 
graphene with massless fermions 15l-ll7l| , this is certainly 
not the case. Some recent experimental works [lil [isl. Il9| 
have demonstrated a relatively high carrier mobility in 



different graphene devices in the high density regime 
(> lO'^cm^V^^s^^) leading a well defined cyclotron reso- 
nance even in the case of small magnetic fields. A clearly 
intriguing problem is to explore how graphene behaves 
when embedded in a cavity resonator. In particular, is it 
possible to achieve ultrastrong coupling between a cavity 
photon and the cyclotron transition in graphene ? If yes, 
is there any qualitative difference with respect to the case 
of massive fermions in semiconductors ? 

In this letter, we present a microscopic theory for the 
cavity QED of graphene under perpendicular magnetic 
field, giving an answer to the questions formulated above. 
Through a quantum field approach (accounting also for 
Coulomb depolarization effects) we show that the ultra- 
strong coupling regime is achievable for graphene. More- 
over, due to the emergence of the Dirac cones with linear 
energy dispersions, the role of the A^ turns out to be 
negligible when the lattice constant is much smaller than 
the photon wavelength. We predict the occurrence of 
a vacuum instability, which is absent in the case of the 
2DEG of massive electrons. 

Graphene is a single layer of carbon atoms forming a 
honeycomb lattice with two equivalent triangular Bra- 
vais sublattices A and B. Each atom of A [B) type 
is connected to its 3 nearest neighbors of B {A) type 
via the displacement vectors ei = ay/2>/2ux + a/2uy, 
62 = —a\/?i/2ux -\- a/2uy and 63 = —aUy. We consider 
here the tight-binding description taking only the first 
nearest neighbors into account with a hopping param- 
eter t. In this case, the first quantization Hamiltonian 
in Q-space can be expressed in a two-components spinor 
basis due to the sublattice isospin[20j. namely: 
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where hr 



iQ = —t X)j=i e~^'^'^^ . In the vicinity of the 
two inequivalent Dirac points (valleys) K and K' marked 
by the two vectors = ±ATT/{3V3a)ux, the low- 

energy excitations are well described by the massless 
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Dirac Hamiltonian T-La.n = (ctKxO'x + f^yf^y) where 
K = Q — K^, ffi {i = x,y,z) are the Pauh matrices, 
vp — {3at)/{2H) is the Fermi velocity lO^m/s) and 
a = ± is the valley isospin index. 

Here, we are interested in the case where a static mag- 
netic field B is applied perpendicularly to the graphene 
plane. This problem can be exactly solved by means 
of the Peierls substitution (which remains valid as long 
as a is much smaller than the magnetic length Iq — 
•y/ {hc)/{eB) in the Dirac Hamiltonian), i.e. replacing 
K with Ho = p + ^Aq where Aq = —ByUx is the vec- 
tor potential (Landau gauge). This yields the graphene 
Landau Levels [Hj whose energies are E^^^ = XhwQ^/n 
where n > 0, A refers to the electron (A = -|-) or hole 
band (A = — ) and loq = vf^/^/Iq- Each Landau level 
has a degeneracy N = QsS/ {2-1:1^), S being the surface 
of the graphene layer and Qs = A accounts for the spin 
and valley isospin degeneracy (we neglect Zeeman and 
valley splittings). Finally we define the filling factor (we 
use a different convention than the usual one[l3, II3|) 

V = pS/N +1/2 (p is the total electron density) so that 

V = n corresponds to the situation where the last fully 
occupied Landau level has an orbital quantum number 
n — 1. In the following, we will study the case of integer 
filling factors and consider for simplicity the case of Fermi 
level in the electron band (A = -I-). Moreover, we will 
consider a graphene layer embedded in a OD cavity res- 
onator. We investigate the coupling between a particular 
cavity photon mode and the transition between the last 
occupied Landau level n = v—1 and the first unoccupied 
one n = V with a frequency wqA;^ = wq (v^ ~ \Jv — l). 
The Dirac fermion field operator in the valley a is written 
as a two-components spinor 
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where the wavefunctions in the valley K are 
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(3) 

The operator c\^n,k,a annihilates a fermion in the state 
labelled by the quantum numbers (A, n, fc, a). C~ = 
\/l — <5n,o/2 and C+ = ^1 + (5n,o/2 are the normaliza- 
tion factors, L is the length of the graphene sheet (taken 
squared for simplicity) along the x and y directions, Lz is 
the cavity length along the z direction and 'Pn,k{y) corre- 
sponds to the eigenfunction of the ID harmonic oscillator 
problem shifted by the guiding centre position yQ — Mq. 
The wavefunctions in the valley K' are obtained by swap- 
ping the spinor ([3]) and changing A — > — A. 

We consider a OD cavity resonator (see Fig. [T]) of vol- 
ume V = LzL"^ that confines the electromagnetic modes 



along the three spatial directions. These cavity modes are 
labelled by the quantized wavevector q = {qx,Clyi<lz) = 



(TTTiaj/L, Trriy/L, Trn^/Lj,), where 



are integer 



numbers. For sake of simplicity, we will consider that the 
cavity length along the z direction is much smaller 
than the cavity transverse size L. This way, we can 
restrict our study to the particular photon mode with 
TLz — 1, neglecting all the higher- lying modes (n^ > 1). 
We will also assume that the graphene layer is placed in 
the middle of the cavity at z = Lz/2 (see Fig. [1]). Of 
course, our theory can be easily generalized to a more 
complicated geometry. The electromagnetic vector po- 
tential reads [221 



/ 27rfc2 



(flq.rjMq,,, + a^,,,Mq,,,) (4) 



where = 1,2 is the photon polarization, lUJcav(g) is the 
cavity frequency, e the cavity dielectric constant and |23| 



Cq^ cos{qxx) sai{qyy) sui{'K z / L z) coaOcoscj) 



Uq.l = 
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yy} oLLLy i\ I ±j z } ^ Sm ( 
—2^/2 s\\Y{qxx) sm{qyy) cos{ttz/ Lz) s\n9 



siTi{qxx) cos{qyy) sm{Trz/Lz) cos 



(5) 



Uq,2 



-Cq^ cos{qxx) sm{qyy) sin(7rz/Lz) sincj) 
Cq^ sm{qxx) cos{qyy) sm{Trz/ Lz) cos0 




(6) 

with cos 6* = qz/\q\ and cos0 — qx/ ^Jqi + 1y- The 
operator aq^n is the photon annihilation operator in the 
cavity mode (q, 77). 

When considering the graphene sheet in the presence 
of a perpendicular static magnetic field embedded in our 
cavity resonator, the electron velocity operator involves 
the total vector potential At — Aq -I- Aem . The coupling 
Hamiltonian is obtained by replacing the momentum p 
in the free Hamiltonian by its gauge invariant form 11 = 
p + I At = IIo -f I Acm- The Hamiltonian ^ becomes 



K^-tY,zf exp 



(7) 



where we have defined the valley-dependent phase fac- 
tor zj"-* = e^*-'^ ■'^J . As we are interested in the quan- 
tum ground state properties and low-energy excitations, 
we can consider the continuum limit (lattice size a — >■ 0, 
which is the relevant limit when a is much smaller than 
any other length scale). In other words, we expand the 
Hamiltonian ([T]) with respect to the lattice parameter a 
and retain the leading contribution. This way, we get 
two terms, the first one is the massless Dirac Hamilto- 
nian 'Hp"'' with static magnetic field, while the second one 



3 




FIG. 1: Left panel: sketch of a cavity resonator of volume V — 
LzL^ (Lz L) embedding a graphene layer with an uniform 
and static magnetic field B perpendicular to the atomic plane. 
Right panel: a fully confined cavity photon mode is supposed to 
be quasi-resonant to the graphene cyclotron transition between 
LLs n — V — 1 and n — u (red horizontal solid lines). 



hI^^ represents the coupling between the Dirac fermions 
and the cavity optical modes: 



early on the bosonic operators Uq^ri and a^^. Conse- 
quently, the total Hamiltonian does not contain any term 
involving the squared vector potential A^^ as it is usu- 
ally the case in cavity QED for massive electrons. Us- 
ing Eqs. ([U, (HI) and the second quantized cou- 
pling Hamiltonian Hi^t — Jd'^r ^f'^ {r)'Hl^^4/{r) becomes 
valley independent so we can omit the index a taking 
the valley degeneracy into account through the factor 
gs. Note that since we are dealing with optical modes, 
we have |qf|/o 1. This condition allows us to neglect 
the LL mixingf23l when considering only the transition 
1/ — > — 1. The final result is 



^^int = ^ fl^q.rj (flq.Tj + Oq,^) (rfq,») + ^q,,,) (10) 

where ilg i = —fig cos 0, rJq,2 = with the vacuum 
Rabi frequency (for ly > I) given by 



(«) 



n 



(a) 



hVF (alio, a; CTj; + no,yO-y) (8) 



Note that the coupling Hamiltonian ^ depends lin- 



afs9s 
47rVe' 



(11) 



The boson annihilation operators dq .fj (polarization 
rj — 1,2) corresponding to the bright collective modes 
are: 



dq,! 



-^^sin L(fc+f )/2 



-^;^cos[,,(fc + |);o^-, 

k 



-^^cos [qy{k^^-^)ll 



cJ,_l JjCi/,/c_gj^ , (12) 

'^i/-l.fe''f,'=-9x • (13) 



The first relevant conclusion of our letter is that 
graphene can enter deeply the ultrastrong coupling 
regime. In fact, in the high filling factor regime (i.e. 
u ^ 1) the cyclotron transition frequency is ljo^d ^ 
ojo/{2y/iy). Hence, the vacuum Rabi frequency normal- 
ized to the transition frequency 



(14) 



can be larger than 1 for large enough filling factors. 
To complete our treatment, we consider the Coulomb 
interaction following a bosonization procedure jisj - Erj l al- 
lowing us to calculate the depolarization shift (Random 
phase approximation (RPA) contribution). The second 
quantized Coulomb Hamiltonian reads: 



(15) 

In the considered geometry, it is convenient to expand 
the Coulomb potential V{r — r') using a 2D Fourier se- 
ries. Finally, we use the same procedure as for the light- 
matter coupling Hamiltonian derivation restricting the 
Coulomb Hamiltonian to describe only scattering pro- 
cesses between LLs v — \ and v thanks to the condition 
\q\l[) ^ 1. This provides the Coulomb Hamiltonian writ- 
ten in terms of the collective modes dq^^ and Cq,,; 



i^coul = ^ flVq-frj {d\ rj + lr,dq,r)^ - {e\ n " 7r,eq,,,)' 



with 7i 



-1, 72 



1 and where Va 



(16) 



32e 



qj_ = {qx,qy) denotes the in-plane wavevector and 
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C^_^\/v — 1 + C^_2V^- Iri the high fihing factor regime, 
we have ~ \f^. 

The Coulomb Hamiltonian reveals the presence of two 
additional dark modes (i.e uncoupled to the electromag- 
netic field) corresponding to two operators Cq^,, and ^. 
Since the bright and dark modes are decoupled by the 
Coulomb interaction, if we look at the ground state prop- 
erties, we can omit the dark modes. After bosonizing the 
kinetic energy |25|, i26j|, the Hamiltonian reads: 



= ^ fu^p{q) m^,^TOq^^ + const. (17) 



where i and ™q,2 corresponds to the niagnetoplas- 
mon modes with energy ujp{q) — y^cuoAi, (wqAi/ + 4Vij). 
Since our bosonic approach allows us to calculate the 
RPA contribution [25|, we find that — >■ when q — >• 0. 
However, our treatment does not include the vertex cor- 
rections (attractive electron- hole Coulomb interaction ) 
nor the difference between the exchange self-energy of the 
electron and that of the hole. For this reason, we do not 
find the renormalization of the transition frequency 
at zero wavevector which is expected since the Kohn's 
theorem[2^ does not apply in graphen eH . Neverthe- 
less, some recent theoretical studies have shown 
that these corrections lead to a relatively moderate renor- 
malization (< 20%) of the Fermi velocity and can thus 
be easily accounted for in our model. In terms of the 
magnetoplasmon modes, the Hamiltonian reads: 



^^tot = hujpiq) ml^rnq^^, + ^cav(g)nl,,^aq,>7 

+ Mq,^ (flq^,, + at ,J (mq^,, + ,^) (18) 



1^0 At, 



where Aqj = ^gCosO^^^ and Aq,2 ^j^. 

In contrast to the Hopfield[3^ quantum Hamil- 
tonian for massive quasi-particles in semiconductor 
microcavitiesjl, such an Hamiltonian is reminiscent 
of the Dicke Hamiltonian, which admits a quantum crit- 
ical point (QCP) beyond which the normal ground state 
becomes unstable. Diagonalizing the Hamiltonian ([18]) 
via the Bogoliubov method, we find that the energy of 
the lower eigenvalue vanishes (existence of a gapless ex- 
citation) when flq - 



^oA„t^cav(q) 



for the branch rj — 1 



and flq — 



4^0 A „ 



for the branch 
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2. Note 



that the no-go theorem for cavity QED quantum phase 
transitions [14| holds as long as the effective kinetic en- 
ergy takes the usual quadratic form /2m* which well 
describes massive electrons in semiconductors (m* is the 
electron band effective mass). In graphene, the peculiar 
crystal structure is essential in determining the electronic 



properties and the effective kinetic energy is governed by 
the massless Dirac Hamiltonian [isi [20| vpp-cr at low en- 
ergy (cr — {(TxT<Jy) being the Pauli matrices). Since this 
Hamiltonian is linear in p, such no-go theorem does not 
necessarily apply and indeed our work (which starts from 
the tight-binding Hamiltonian of graphene) predicts that 
a vacuum instability occurs. 
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FIG. 2: Left panel: normalized frequencies aj-,,,,/ajcav (q) 
(solid lines) of the graphene magnetopolariton lower and 
upper branches as a function of the carrier density for a 
given magnetic field B = 50 mT. The dashed lines depict 
the bare cavity and magnetoplasmon frequencies. Parame- 
ters: e = 3.9, Lz = 760 /^m and L = 5Lz, optical mode 
with q = {ux = l,ny — 1,71^ = 1) and ujq — 0.628 THz rad~^ 
(/q = 0.1 THz). Right panel: same quantities (unnormalized 
and divided by 2-k) as a function of the magnetic field just 
below the critical density pc = 1.4 x lO^'^cm"^. 

The left panel of the Fig. [2] shows the frequencies 
of the magnetopolariton excitations ijJj^n, (j = LP, UP) 
normalized to the cavity mode frequency versus elec- 
tron density. We see that by increasing the density, 
the lower branch rj — \ has a nonmonotonical behav- 
ior and vanishes when reaching a critical value pc — 
1.4x lO^^cm^^. On the right panel, the magnetopolariton 
frequencies — a;j_,,/27r are depicted as a function of 
the magnetic field just below the critical density pc- We 
notice a strongly asymmetric dispersion when increasing 
the magnetic field together with the vanishing frequency 
of the lower branch /lp,i. These characteristics are a sig- 
nature of such a vacuum instability. In order to describe 
the properties of the Dicke Hamiltonian above the QCP, 
an Holstein-Primakoffjs^l bosonization procedurefssjl can 
be applied. However, such a mapping does not hold in 
this case because of the non- vertical excitations at finite 
wavevector [2^ , making the analytical description of the 
excitations above the QCP highly non-trivial and will be 
the scope of future investigations. 

In conclusion, we have shown that the peculiar crystal 
structure of graphene leads to remarkable cavity QED 
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properties. In particular, the coupling of the graphene 
cyclotron transition is qualitatively different for the case 
of massive electrons in semiconductors. In the ultra- 
strong coupling regime, a vacuum instability analogous 
to the one occuring in the Dicke-model can also occur 
for graphene, whereas for massive quasi-particles this is 
not the case. In free space graphene (without cavity), 
the anomalous properties of Dirac quasiparticles have 
been shown to lead to an unusual behavior of the op- 
tical conductivities in the presence of a perpendicular 
magnetic field [sij. For cavity embedded Dirac fermions, 
our work paves the way to interesting developments, such 
as cavity-controlled magnetotransport in graphene. We 
would like to thank M. O. Goerbig, S. De Liberato and 
P. Nataf for fruitful discussions. 
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